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The QCD pressure is a most fundamental quantity, for which lattice data is available 
up to a few times the critical temperature T c . Perturbation theory, even at very high 
temperatures, has serious convergence problems. Combining analytical and 3d numerical 
methods, we show that it is possible to compute the QCD pressure from about 2T C to 
infinity. We also describe an algorithm to generate and classify high order Feynman 
diagrams which is tailored to minimize computational effort. 



Introduction. The properties of QCD matter are expected to change above a critical 
temperature of the order of 200 MeV. While the low-temperature phase is governed by 
bound states, such as mesons, the high-temperature phase should, due to asymptotic 
freedom, look more like a gas of free quarks and gluons. Any observable witnessing this 
change is therefore a potential candidate for (direct or indirect) measurements in heavy- 
ion collision experiments. 

One such observable clearly is the free energy density of the system. The rough picture is 
that it is, according to the Stefan-Boltzmann law, proportional to the number of effective 
degrees of freedom. For vanishing baryon density ^ = and at temperature T, the free 
energy density of QCD is simply given by the functional integral 



/ = -Iln|p[A^/]exp^-^ 1/T dr|^ 3 x 



(i) 



Note that, in the thermodynamic limit of infinite volume V, the pressure p of the plasma 
is given directly by p = —f. Below, we will choose to display results for the pressure. 

The most direct way to evaluate this integral would now be to measure it numerically 
on the lattice. In fact, this has been done by a number of groups. While the results for 
Nf = are rather complete, they are rapidly developing for a finite number of fermion 
flavours Nf (e.g. . The general picture emerging from these lattice simulations is the 
following: Normalizing the pressure to zero below T c (in practice, one can only measure 
derivatives of the free energy, leaving an integration constant to be fixed), it rises sharply 
in the interval (1 — 2)T C , to level off at a few times T c . At the highest temperatures used in 
the simulations, typically (4 — 5)T C , the deviation to the Stefan-Boltzmann limit is of the 
order of 15%. This general picture is surprisingly stable with respect to different values 
of Nf. At even higher temperatures, the pressure is then expected to asymptotically 
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approach the ideal-gas limit p (T) = (tt 2 T 4 /45)(A^ c 2 -1 + (7/4)N c N f ), where N c denotes 
the number of colours. 

It turns out that this deviation of 15% is too large to be understood in terms of ordinary 
perturbation theory. In a series of impressive works, the expansion has been driven to 5th 
order in the gauge coupling g 0, 

221 = l + c 2 g 2 + c 3 g 3 +(c',\ng + Ci )g 4 + c 5 g 5 + 0(g 6 \ng,g 6 ). (2) 
) 

While the coefficients are known analytically, convergence properties are extremely poor, 
cf. Fig. [I], at least for all physically relevant temperatures. 

Facing the poor convergence of the perturbative series, in the past few years a lot of 
effort has gone into refined and/or alternative approaches, in order to gain an analytic 
understanding of the high-temperature behaviour of the pressure. The spectrum ranging 
from Pade-Borel resummations over using effective masses to employing selfconsistent 
approximations ||, a general feature of these works is the suppression of infrared (long- 
distance) effects. While this suppression does not seem to be crucial in the computation 
of the pressure (which appears to be a short- distance dominated observable)Q, the aim 
of this talk is to review our framework to resum the long-distance contributions to the 
pressure to all orders ||. 

Combined analytic and 3d numerical method. A way to understand the poor con- 
vergence of the ordinary perturbative expansion is the observation that at small gauge cou- 
pling g, the system undergoes dimensional reduction (see e.g. |J and references therein). 
The scale hierarchy gT <C ttT allows to perturbatively construct an effective theory for 
the "soft" modes (momenta oc gT) by integrating out the "hard" modes (oc ttT). In the 
case of QCD, this effective theory is a 3d SU(iV c ) + adjoint Higgs model: 

Ad = ^TF^ + ^i^A^ + ^mlTTAl + ^XA^iAlf + dC^. (3) 

While the last term represents higher-order operators, which we do not take into account 
here because their contributions are parametrically suppressed |7|], the parameters of the 
first terms (g 3 , m 2 Dl \a) are related to the physical parameters of the full 4d theory (T, 
Ajyfg). Using optimized next-to-leading order perturbation theory ]J[ (let us introduce 
dimensionless parameters^ x, y and set Nf = 0, N c = 3 here), they read: 

g\ _ 8tt 2 /11 _ \ A _ 3/11 _m 2 D 3 9 

T ~ ln(6.742T/A M g) ' X ~ ~gj ~ \z(5.371T/A m ) ' V ~ It ~ 8^ + 16^ ' U 

Two comments are now in order. First, can be used to reproduce Eq. (0) in a 
technically simple way ||. This shows that the bad convergence is due precisely to the 
soft degrees of freedom, and it provides a clearer understanding of the contributions from 
separate physical scales. 

Second, the effective theory is confining, hence non-perturbative |§. This fact directly 
leads us to the conclusion that the only way to systematically include the long-distance 
contributions to the pressure is to treat on the lattice. 

fact, the last method mentioned connects to the lattice data available quite nicely from above. 
2 The 3d gauge coupling g\ has the dimension of a mass. 
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To proceed, we rewrite the pressure, up to hard-scale g 6 contributions, as 



P(T) _ 5x _ 45_ (g\ 

Po{T) 2 8tt 2 \T 



(5) 



where 5 ~ 10~ 4 and the dependence on the scale p,^, which originates from an infrared 
divergence of the 4d part, cancels against a similar ultraviolet termP] in the dimensionless 
3d free energy density jFj^g(x, y), 

which should be measured on the lattice. This requires a measurement of the quadratic 
as well as quartic adjoint Higgs field condensates (which determine the partial derivatives 
under the integral), as well as a (4- loop) perturbative computation in lattice regulariza- 
tion, to match to the MS scheme. We choose to fix the integration constant perturbatively 
at high temperatures T ~ 10 n T c , where one is confident that the expansion converges^. 

Before displaying first results for the pressure obtained by the above strategy, let us 
illuminate one aspect in more detail here, namely the setup of its high-order perturbative 
expansion in an efficient manner. 

Diagrammatics. The enumeration of Feynman diagrams contributing to a specific loop 
order together with a derivation of the accompanying symmetry factors seems to be the 
'trivial' part of any perturbative calculation. At higher orders, it is our experience that an 
efficient algorithmic setup of this initial step proves necessary. This does not only assure 
completeness of the required set of diagrams, but it also has the potential of streamlining 
the subsequent integration step considerably by grouping together sets of related diagrams, 
thereby avoiding an unnecessary repetition of subdiagram computations ||. 

The main idea is to utilize the very efficient notion of skeleton (2-particle-irreducible, 
2PI) diagrams to achieve the above-mentioned grouping (into skeleton and ring diagrams, 
see below). For simplicity, we will turn to a generic bosonic 3 + 4 theory here. The 
skeleton expansion for the free energy as a functional of the full propagator D reads [|l(| 



- F[D] = -- (Tr In IT 1 + Tr U[D]d) + $[D] . (7) 

Here $[D] collects all 2PI vacuum diagrams. The full propagators D are related to their 
corresponding self-energies by D" 1 = A -1 — II where A are the free propagators. The 
partition function has an extremal property, such that the variation of F with respect 
to the full propagator vanishes, giving a relation between skeletons and self-energies, 
8dQ[D] = | II[Z)]. Pictorially, this corresponds to obtaining a self-energy by "cutting 
a propagator" in all possible ways in the set of vacuum skeletons. Hence, knowing the 
skeletons alone provides full information. 

To utilize Eq. (0), we need to rewrite it as a loop expansion, which can be achieved 
by a straightforward Taylor expansion around the free propagator A. The result, up to 



3 This is precisely the way the effective theory is set up: dependence on a matching scale has to cancel. 
4 T c is of the order of A^g, the coefficient can be measured on the lattice. 
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4-loop order, reads diagrammatically (denoting by Fo the non-interacting result) 
_jr = _jr o + ^[A] + (* 3 [A]+i<0>) + ^ , (8) 

clarifying our terminology of having achieved a separation into skeleton and ring diagrams. 
The notation here is that 3> n [A] are n-loop skeletons built up of free propagators, while 
a circle/square with n inside denotes what we term an n-loop irreducible/reducible self- 
energy n" r /II„ . More precisely, the irreducible self-energies derive directly from the 
skeletons, while reducible self-energies have (at leastQ) one further self-energy insertion, 

-©- =n-[A] = 2d A $ n+1 [A] , =n^[A] = njr_ 1 [A + An«A]| .(9) 

n— loop 

Hence it is now clear that Eq. @ expresses the loop expansion of the free energy in an 
economic way in terms of the irreducible skeletons $ n [A]: either as direct contributions, 
or as self-energy insertions, obtained from the same skeletons by derivatives. 

While everything is reduced to knowing the skeletons and their symmetry factors now, 
the key advance is that one can write down a closed exact equation which directly gener- 
ates the n-loop skeletons 0: 

*«i A '-;rM^O + i(X) + KD^(©L- n> - 2 - (10) 

Lines are bare propagators A, and vertices are bare and full ones (shaded blob), respec- 
tively. The full vertices are generated in the standard way by the tower of Schwinger-Dyson 
equations, which we do not repeat here. For a functional derivation as well as a graphical 



representation see, e.g., 11 



It is now straightforward to utilize the concept of Eqs. (|8"D , ([TD|) to generate all diagrams 
needed for the theory in Eq. (|3|), to serve as a starting point for the (up to) 4-loop 
computations needed. 

Discussion. As a first result obtained along the strategy outlined above, Fig. ^| shows the 
normalized pressure. The integration constant has been fixed perturbatively on the 3-loop 
level, allowing for an additional constant en., which represents an (up to now) unknown g® 
contribution. In principle, this constant can be determined in a setup equivalent to the 
above, after splitting off its perturbative part: A further reduction step relates en to the 
free energy of 3d pure gauge theory, which is at the core of the famous non-perturbative 
g 6 term, but can nevertheless be determined on the lattice. On the lattice side, we have 
only included the quadratic scalar condensate, while at temperatures closer to T c , the 
quartic one will become important as well. 

While more work is required (and in progress), we wish to point out two important 
trends seen in Figs. First, the outcome is sensitive to the non-perturbative parameter 
en, which in principle can be determined by additional computations. Clearly, there exists 
a range for that parameter (say, O(10)), which leads to a sensible result. 



5 In higher orders, reducible self-energies divide naturally into classes according to the number and type 
of self-energy insertions. At 4-loop level there is only one class, so we do not introduce extra indices here. 
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Figure 1. Perturbative pressure vs. lattice 
data for Nj = (from Q). The notation 
is like in Eq. (|), corresponding to p up to 
g 2 , g 3 , g 4 and g 5 . 



Figure 2. The pressure after inclusion of 
the long-distance part according to Eq. @ 
(from Statistical errors are shown 

only for e = 10. 



Second, at T > 30A^g the curves for g 5 (yi ) and eo = 10 fall almost on top of 
each other, signalling a cancellation of all higher-order terms (determined here by the 
quadratic Higgs condensate) against the large non-perturbative g 6 contribution. Hence, 
in this temperature range the pressure is indeed dominated by short- distance effects. 

While we have mostly presented results for QCD at Nf = and zero baryon number, 
inclusion of Nf fermion flavours as well as a baryon chemical potential //& pose no further 
complications, and hence provide for a natural extension of this investigation. 
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